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Abstract 

We obtain series expansion formulas for the Hadamard fractional integral and fractional 
derivative of a smooth function. When considering finite sums only, an upper bound for the 
^f) ' error is given. Numerical simulations show the efficiency of the approximation method. 
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1 Introduction 

In general terms, Fractional Calculus allows to define integrals and derivatives of arbitrary real or 
complex order, and can be seen as a generalization of ordinary calculus. A fractional derivative 
of order a > 0, when a is integer, coincides with the classical derivative of order n 6 N, while a 
fractional integral is an n-folded integral. Although there exist in the literature a large number 
of definitions for fractional operators (integrals and derivatives), the Riemann-Liouville and 
Caputo are the most common for fractional derivatives, and for fractional integrals the usual 
one is the Riemann-Liouville definition. In our work we consider the Hadamard fractional 
integral and fractional derivative. Although the definitions go back to the works of Hadamard 
in 1892 [5], this type of operators are not yet well studied and much exists to be done. 
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As well-known, for most problems involving fractional operators, such as fractional differen- 
tial equations or fractional control problems, one cannot provide methods to compute the exact 
solutions analytically. Therefore, numerical methods need to be employed. Typically, problems 
are addressed by looking to the fractional operators as special types of integrals and, using par- 
titions of the domains, writing them as finite sums, with some error for the final result. Our 
approach is distinct from this, in the sense that we seek expansion formulas for the Hadamard 
fractional operators with integer-order derivatives. In this way we can rewrite the original prob- 
lem, which depends on fractional operators, as a new one that involves integer derivatives only 
and then successfully apply the standard methods to obtain the desired solution. Also, in some 
cases, it can be easier to find the fractional derivative or integral using those expansions, instead 
of applying the direct definitions. In contrast with [2], where an expansion for the Riemann- 
Liouville fractional derivative is given, our expansions for the Hadamard fractional derivative 
and integral do not omit the first derivative, allowing one to obtain considerable better accuracy 
in computation. 

The paper is organized in the following way. In Section [2] we review some concepts of 
fractional calculus. Decomposition formulas for the left and right Hadamard fractional integrals 
are given in Section [31 together with approximation formulas and error estimations. Following 
the same approach, similar formulas are obtained for the left and right Hadamard fractional 
derivatives in Section 2J In Section [5] we test the efficiency of such approximations with some 
examples, comparing the analytical/exact solution with the numerical approximation. 

2 Preliminaries 

In this section we review some necessary definitions for our present work, namely the Hadamard 
fractional integral and derivative. For more on fractional calculus, we refer the interested reader 
to [U [101 QTJ [13]. For related work on Hadamard fractional operators, see [21 El El l9l [12]. 

Let a, b be two reals with < a < b and x : [a, b] — > M be an integrable function. The left 
and right Hadamard fractional integrals of order a > are defined by 



respectively. These integrals were introduced by Hadamard in [5] in the special case a = 0. 
When a = m is an integer, these fractional integrals are m-folded integrals (see, e.g., [31 [6]): 




and 
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and 

«. xW = /**»/•**.../• fW*., 

For fractional derivatives, we also consider left and right operators. For a > 0, the left and 
right Hadamard fractional derivatives of order a are defined by 

a V?x(t) = (41 " f f In 1 ) ^ —dr, t e]a, b[, 



dt J r(n — a) J a \ t 
and 

t VZx(t)= (-t±Y-J— T /Vln^^^r, te] 0> 6L 



/ r(n — a) J t \ t 
respectively, with n = [a] + 1. When a = m is an integer, we have (cf. |5]) 

(j \ 771 / rl \ m 

t-j x{t)^d t V^x{t)=[-t-\ x(t). 

Hadamard's fractional integrals and derivatives can be seen as inverse operations of each 
other (see Property 2.28 and Theorem 2.3 of [8]). 

When a £ (0, 1) and x 6 v!C[a, 6], with ^4C[a, 6] denoting the space of all absolutely contin- 
uous functions x : [a, 6] — >■ R, the Hadamard fractional derivatives may be expressed by 



and 



For an arbitrary a > we refer the reader to Theorem 3.2]. If a function x admits derivatives 
of any order, then expansion formulas for the Hadamard fractional integrals and derivatives of 
x, in terms of its integer-order derivatives, are given in [H Theorem 17]: 

oo 

l?x(t) = Y,S{-a,k)t k x ( - k \t) 

k=0 

and 

oo 

oV?x(t) = Y,S(a,k)t k x^(t), 

k=0 

where 

k 



is the Stirling function. 
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3 An expansion formula for the Hadamard fractional integral 

In this section we consider the class of differentiable functions up to order n + 1, x £ C n+1 [a, b], 
and deduce expansion formulas for the Hadamard fractional integrals in terms of x^(-), for 
i G {0, . . . , n}. Before presenting the result in its full extension, we briefly explain the techniques 
involved for the particular case n = 2. To that purpose, let x <E C 3 [a,b]. Integrating by parts 
three times, we obtain 



1 



1 



, In - 

r(a) J a T\ T 
1 f. t 



t 



a-1 



r(a + i) 
1 



In 



x{t)cLt 

1 /"* 1 / t\ a 
x(a) — — r / — (in— rxfrW? 



In - I x(a) + 
a 



1 



r(a + i) 

-f^)/'-K ln 7)° +, ^w + 



a+1 



1 



r(a + i) 
1 



r(a + 2) V a 

a+l 

T 2 x(j))dT 
a+l 



ax (a) 



In 



r(a + 3) V « 



i \ 1 

In - I x(a) + — > i ■ ■ ■ 

a J v ' r(a + 2) V a 

A a+2 

In- I (ax(a) + a 2 x(a)) 



ax(a) 



+ 



1 



In 



a+2 



T(a + 3) J a \ T 

On the other hand, using the binomial formula, we have 



(x(t) + 3tx(t) + t 2 x(t))(It 



In 



a+2 



= In 



In 



a+2 



1 - 



In- 



t \ a+2 



In 4 



/r+ 2 ^ I>-a- 2) (In I)' 



2)p! (In if 

This series converges since r £ [a, i] and a + 2 > 0. Combining these formulas, we get 

a+l 



1 



+ 



r(a + 3) 



r(a + l) 



In* 

a 



In-) s(a)+ 1 (ln- 
a / r(a + 2) \ a 

a+2 oo 



ax(a)+- 



1 



,a+2 



V- 



r(p - a - 2; 



r(a + 3) \ a 
-p-p I (in -V (x{t) + 3rx(r) + r 2 x (r)) dr 



(ax(a)+a x(a)) 



p= - v -a - 2)pl (ln|) 

Now, split the series into the two cases p = and p = 1 ... 00, and integrate by parts the second 
one. We obtain 



a T?x{t) 



1 



r(a + l) 

1 



r(a + 3) V a 



t\ 1 

In - x(a) + — r 

a J w T(a + 2) 

,x a+2 

ln M nun , r> 



In* 

a 



a+l 



(tx(t) + t 2 x(t)) 



In* 



a+2 00 



ax(a) 

r(p - a - 2) 



r(a + 2) V ay ^r(-a-l)(p 



^T(-a-2)p! 
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Repeating this procedure two more times, we obtain the following: 



1 



r(a + l) V a 



In- x(t) 



p=3 



r(p - a - 2) 
r(-a)(p-2)! 



+ 



1 



■ In — 

r(a + 2) V a 



a+l 



p=2 



r(p - a - 2) 
r(-a-l)(p-l)! 



1 



In 



Q + 2 



T(a + 3) V a 



(fcc(£) + t 2 x{t)) 



r(p - a - 2) 



r(p - a - 2) 
, T(-a - 2)pl 

1 T\P-3 X (t) 



r(a)V ^r(-a + l)(p-3)!(lni) p 7 a 

or, in a more concise way, 



In 



dr, 



aTfxit) = A (a) In - x(t) + Ai(a) In 



Q+2 



a+l 



tx(t) 



+A 2 (a) (in- j (ii(t) + t 2 i(t)) + ^ B(a,p) ( In - 

p=3 



a+2-p 



V p (t) 



with 



A (a) 
A x [a) 
Ai(a) 



1 



r(a + i) 
1 

" r(a + 2) 
1 

" r(a + 3) 
B(a,p) 



oo 



p=2 



i+E 



r(p - a - 2) 
r(-a)(p-2)! 

r(p - a - 2) 
r(-a-l)(p-l)! 



r(p - a - 2) 



and 



V P (t)= [ (p-2)(lnl 



^T(-a-2)p! 

r(p - a - 2) 
r(a)r(l -a)(p-2)! 

p-3 ifrl 



-dr, 



where we assume the series and the integral V p to be convergent. 



(3) 
(4) 



Remark 3.1. When useful, namely on fractional differential equations problems, we can define 
V p as in @ by the the solution of the system 



V p (t) = (p-2) (lni) 
V p (a) = 



t \P~3 z(t) 



/or aZZ p = 3, 4, . . . 
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We now discuss the convergence of the series involved in the definitions of Ai(a), for i 6 
{0, 1, 2}. Simply observe that 

T(p-a- 2) 



E 



T{-a-i){p-2 + i)\ 



p=3— i 

and 1^0(0,2;) converges absolutely when \x\ = 1 if a < ([1, Theorem 2.1.2]). 

For numerical purposes, only finite sums are considered, and thus the Hadamard left frac- 
tional integral is approximated by the decomposition 

o+l 



al?x{t) ^A (a,N) In- x(t) + At{a, N) In 



t 

a J 



tx(t) 



.\Q+2 N , sa+2-p (5) 

-A 2 (a,N)(\n-J (tx(t) + t 2 x(t)) + Y^B {a, p) (in- j V p {t) 



with 



A (a,N) 
Ai(a,N) 
A 2 (a,N) 



1 



r(a + i) 
1 

T(a + 2) 
1 

T(a + 3) 



1 | y r(p-a-2) 
^ 3 r(-a)(p-2)! 

i + ^ r(p-«-2) 

p=2 

AT 



r(-a-l)(p-l)! 
r(p - a - 2) 



^ r(-a - 2) P ! 

B(a,p) and V^(t) as in (J3J) — (J31) , and iV > 3. We proceed with an estimation for the error on 
such approximation. We have proven before that 



a l?x(t) 



1 



• In — 

r(a + l) V a 



x(a)+ 



r(a + 2) 



In- 



0+1 



ax(a)+ 



1 ^ n 

\ 1 ■■ + 3) \ a 



a+2 



(ax(a)+a 2 x(a)) 



1 



ln- 



0+2 r t 00 



E 



r(p-a-2) Qnj 



r(a + 3) V a/ 7a ^T(-a-2)p! 
When we consider finite sums up to order iV, the error is given by 



x(t) + 3tx(t) + r 'x(t))(It 



\Etr(t)\ 



1 



with 



r(a + 3) 
Rn{t 



In 



a+2 /-t 



E 

p=iV+ 



/ R n (t)(x(t) + 3tx(t) + t 2 x{t))cLt 

j a 

r(p-o-2) (ln^) P 



t \p ' 



^(-a-2)^! (lni) 



Since r £ we have 



iflw(T)i < £ 



a + 2 
P 



s E 

p=7V+l 



£2. e (a+2) 2 +a+2 



,a+3 



p=Af+l 

/■oo e ( a +2) 2 +a+2 g (a+2) 2 +a+2 

"iv ^ = (a + 2)Ara+2- 



6 



Therefore, 

\E tr (t)\ < 
where 



/ f\ a+2 e ("+2) 2 +a+2 

V ln a) (a + 2)N^ ^ " a)Ll{t) + 3{t ~ + {t ~ a ^ L ^ ' 



r(a + 3) V a J ( 



Li(t) = max |x w (r)|, i G {1,2,3}. 

re[o,t] 



Following similar arguments as done for n = 2, we can prove the general case with an 
expansion up to the derivative of order n. First, we introduce a notation. Given k € N U {0}, 
we define the sequences x kj o(t) and x kt i(t) recursively by the formulas 

XQfl{t) = x(t) and x k+ i,o(t) = t^x kj0 (t), for k G N U {0}, 



and 



Xo,i{t) = x(t) and x fc +i,i(t) = ^(tx k>1 (t)), for k £ N U {0}. 



Theorem 3.2. Let n £ N, < a < i and x : [a, b] — > R 6e a function of class C n+1 . Then, 

a+i oo / ^ \ a+n— p 

a(a,p) i m ' 

p=n+l 



a Xfx(t) = ^^(a) (in 



i=0 



: ,o(t)+ 2 S(a,p)^Iii-J F p (t) 



with 



Ai(a) = 
B(a,p) 



r(a + i + 1 

F(p — a — n 



» + E 



r(p — a — n) 



' r(— a — i)(p — n + i)l 



r(a)r(l - a)(p-n)V 
£(p-n)(lnl)—^) dr . 



Moreover, if we consider the approximation 



a+i 



N 



a Zfx(t)^^(a,iV) (lni) z i)0 (t) + £ 5(a,p) f In ^ 

p=n+l 



a+n—p 



v P (t) 



with N > n + 1 and 

A (a, AT) 



1 



v 



T(a + i + 1 
f/ien i/ie error is bounded by the expression 

\E tr (t)\ < L n (t) 

where 



i+ E 



V{p- a- n) 



' T(—a — i)(p — n + i)\ 



e (a+n) 2 +a+n / £ 

r(a + n + l)(a + n)iV Q+n V 11 a 



a+n 



(t-a), 



L n (t) = max |x„,i(t)|. 

r£[a,t] 
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Proof. Applying integration by parts repeatedly and the binomial formula, we arrive to 

a+i 



a T?x{t) =J2 



i=0 



T(a + i + 1 
1 



In* 

a 



Xi,o{a) 



^\«+n oo^ T(p — a — n 
T{a + n + l)\^a) ^ T(-a - n)p\ (lni) 



To achieve the expansion formula, we repeat the same procedure as for the case n = 2: we split 
the sum into two parts (the first term plus the remainings) and integrate by parts the second 
one. The convergence of the series Ai{a) is ensured by the relation 



r(p — a — n) 



FY— a — i)(p — n + i)\ 

p=n—i+l 



The error on the approximation is given by 

1 



Fo(-a-iA) - 1. 



l-TQ 



\Etr{t)\ = 



T(a + n + 1) 



In- 



a+n 



\ R N {r)x n ±{T)dn 

J a 



with 



Also, for r 6 [a, t], 



\Rn(t)\< £ 

p=N+l ^ 



; v+i r(-a-n)p! (In if 

(a+n) 2 +a+n 



a + n 



< 



(a + n)iV' 



a+n 



□ 



We remark that the error formula tends to zero as TV increases. Similarly to what was done 
with the left fractional integral, we can also expand the right Hadamard fractional integral. 

Theorem 3.3. Let n £ N, < a < i and x : [a, b] — > R be a function of class C n+1 . Then, 
t lZx(t) = ^Ma)(ln-j x ifi (t)+ B(a,p)hn-\ W p (t) 

V / p=n+l ^ ' 



with 



i=0 



Ai{a) 

B(a,p) 
W p (t) 



(-iy 



T(p - a - n) 



V{a + i + \) 
T(p — a — n 

T(a)r(l - a){p-n)V 
rb / fo\p-n-l 



' T(—a — i)(p — n + i)\ 

p=n—t+l 



n) ( In — 



X(T 



dr. 



8 



Remark 3.4. Analogously to what was done for the left fractional integral, one can consider 
an approximation for the right fractional integral by considering finite sums in the expansion 
obtained in Theorem \3.3l In this case, the error is bounded by 



\E tr (t)\ < L n {t) 



where 



e (a+n) 2 +a+n / ^ 

r(a + n + l)(a + n)N a+n \ U t 



L n (t) = max |x„,i(t)|. 

re[t,b] 



a+n 



(b-t), 



4 An expansion formula for the Hadamard fractional derivative 

Starting with formulas (P) and ([2]) , and applying similar techniques as presented in Section [3l 
we are able to present expansion formulas, and respective approximation formulas with an error 
estimation, for the left and right Hadamard fractional derivatives. Due to restrictions on the 
number of pages, we omit the details here and just exhibit the results. 
Given n G N and x G C n+1 [a, b], we have 



a Vfx(t) 



r(i -a) 

oo 

+ E 

p=n+l 



t\~" - ft 

In- ] x(t) + J2Ai(a) (in 

i=l 



\1 - a)T(a)(p - n)\ 



in* 

a 



x(t) 



with 



Ai(a) 



B(a,p) 



T(i + 1 - a) 



i+ E 



T(p + a — n) 



p=n—i+ 



T(p + a — n) 
T(a — i){p — n + i)\ 



ie {l,...,n}, 



-, pe{n + l,...}, 



T(-a)T(l + a)(p-n)\ : 
V P it) =r a (p-n)(\n^) p - n - 1 ^dr, P e{n + 1,...}. 



When we consider finite sums, 



i=0 



v 



a V?x(t)*sJ2M<x>N) (in-) x i)0 (t)+ B(a,p)[bx- 

p=n+l 



n — ol — p 



v P (t) 



with 



Ma,N) 



1 



r(» + i-a) 

and the error is bounded by 

\Etr(t)\ < L n (t) 



N 



» + E 



p=n—i+l 



r(p + a — n) 
T(a — — n + i)! 



i G {0,...,n}, 



3 (n— a) 2 +n— o 



r(n + l-a)(n-a)iV n - 



In* 



(t-a), 
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where 



L n (t) = max |x n i(r)|. 

re[a,t] 



Remark 4.1. The series involved in the definition of Ai are convergent, for all i £ {1, . . . , n}. 
This is due to the fact that 



E 



Y{p + a — n) 
r(a — i)(p — n + i)! 



iF (a - i, 1) - 1 



p=n— i+l 

and iFo(q — i, 1) converges (If, Theorem 2.1.1]), since a — i < 0. 

Remark 4.2. For i/ie rig/i£ Hadamard fractional derivative, the expansion reads as 



1 



• In 

I - a) V t 



+ E 

p=n+l 



B(a,p) ( ln~ 



x{t) + ^Ai(a) ( In 

i=l 

i \ n— a— p 

W p (t) 



^i,o(*) 
T(p + a — n) 



r(l -a)r(a)O-n)! V * 



x(t) 



with 



-1) 



T(i + 1 - a) 



+ E f 



T(p + a — n) 



p=n— i+1 



T(p + a — n) 
(a — i)(p — n + i)! 



, ie{l,...,n}, 



r(-a)r(l + a)(p-n)! 

6 / l\P-i— 1 

(p — n) ( In 



, pG{n+l, ...}, 



A(a) = 

S(a,p) = 
W p (t) = 

5 Examples 



We obtained approximation formulas for the Hadamard fractional integrals and derivatives, and 
an upper bound for the error on such decompositions. In this section we study several cases, 
comparing the solution with the approximations. To gather more information on the accuracy, 
we evaluate the error using the distance 



dist 



(X(t) -X{t)) 2 dt, 



where X(t) is the exact formula and X(t) the approximation. To begin with, we consider a = 0.5 
and functions x\{t) = hit and X2(t) = 1 with t £ [1, 10]. Then, 



/ — 3 — 

lX *°' 5xi(t) = f^sy and iX °' 5x2 ^ 



Vint 

rp) 
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(cf. [H Property 2.24]). We consider the expansion formula for n = 2 as in © for both cases. 
We obtain then the approximations 



i^xi (t) 



N 



A (0.5, N) + Ai(0.5, AO + V 5(0.5,p)^4 



V\n 3 



and 



iX t a5 x 2 (t) 



jV 



A (0.5,N) + Y, 5(0.5, p) 

p=3 



Oni. 



The results are exemplified in Figures 1(a) and 1(b) As can be seen, the value A^ = 3 is enough 



in order to obtain a good accuracy in the sense of the error function. 





(a) i2?- 5 (lnt) 



(b) !2?' 5 (1) 



Figure 1: Analytic vs. numerical approximation for n = 2. 



We now test the approximation on the power functions x^{t) = t 4 and Xi(t) = t 9 , with 
t e [1,2]. Observe first that 

-0.5 

| 7-O.r,, ,/,• , 1 / / ,. 1 \ 



In- 



r(0.5) 

by the change of variables £ = In -. In our cases, 
0.8862269255 



r k - l dr 



t K 



r(o.5) J 



in* 



r(o.5) 



4 p, a — s , f-rO ^ / Qn 0.5908179503 q 
t 4 erf(2Vlni) and i2?' 5 (i 9 ) « 



r(o.5) 



rerf(3Vlnt), 



where erf(-) is the error function. In Figures 2(a) and 2(b) we show approximations for several 
values of N. We mention that, as A" increases, the error decreases and thus we obtain a better 
approximation . 

Another way to obtain different expansion formulas is to vary n. To exemplify, we choose the 
previous test functions Xi, for i = 1,2, 3, 4, and consider the cases n = 2, 3, 4 with A^ = 5 fixed. 
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(a) il t a5 (i 4 ) (b) 1 l? 5 (t 9 ) 



Figure 2: Analytic vs. numerical approximation for n = 2. 



The results are shown in Fig ures |3(a)[ |3(b)[ |3(c)| and |3(d)[ Observe that as n increases, the error 
may increase. This can be easily explained by analysis of the error formula, and the values of 
the sequence £(fc 5 o) involved. For example, for X4 we have X(k,o){t) = 9 k t 9 , for k = . . . ,n. This 
suggests that, when we increase the value of n and the function grows fast, in order to obtain a 
better accuracy on the method, the value of N should also increase. 

We now proceed with some examples for the Hadamard fractional derivatives. The test 
functions are the same as before, and in Figures 4(a), |4(b) 4(c) and |4(dJ| we exemplify the 
results. In this case, 



iV?- 5 x 4 (t) 



Vint 



r(o.5jvhTt' 
r(o.5|vln7 



r(o.5)vint 



+ 



0.8862269255 

rpj 

0.5908179503 



T(0.5) 



4t 4 erf(2Vlnl), 
9i 9 erf(3Vlnl). 



One main advantage of this method is that we can replace fractional integrals and fractional 
derivatives as a sum of integer /classical derivatives, and by doing this we are rewriting the 
original problem, that falls in the theory of fractional calculus, into a new one where we can 
apply the already known techniques (analytical or numerical) and thus solving it. For example, 
when in presence of a fractional integral or a fractional derivative, with a number of initial 
conditions, replace the fractional operator by the appropriate approximation, with the value of 
n given by the number of initial conditions. 
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Figure 3: Analytic vs. numerical approximation for n = 2, 3, 4 and N = 5. 



For example, consider the problem 



x(l) = 0. 

Obviously, x(t) = hit is a solution for ([6]). Since we have only one initial condition, we replace 
the operator \D®^(-) by the expansion with n = 1 and thus obtaining 



N 



[1 + A, (0.5, TV) (lnt)- - 5 ] x(t) +A l (0.5,N)(lnt)°-hx{t) + ^ 5(0.5,p)(lnt) a5 - p Vp(t) 



V p (t) = (p- l)(lntr 2 ^, p = 2, 3, ... ,N, 
x(l)=0, 

V P (1) = 0, p = 2,3,...,N. 



p=2 



r(i.B) 



+ lnt, 



(7) 
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(a) i2? t °- 8 (lnt) 



(b) !© t - 5 (l) 




1 1.2 1.4 1.6 1.8 2 1 1.2 1.4 1.6 1.8 2 

t t 



(C) lP?' 5 (t 4 ) (d) !V ( l 5 (t 9 ) 

Figure 4: Analytic vs. numerical approximation for n = 2. 

In Figure [5] we compare the analytical solution of the FDE ([6]) with the numerical result for 
N = 2 in (0). 
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Figure 5: Analytic vs. numerical approximation for the FDE © with one initial condition. 

References 

[1] G. E. Andrews, R. Askey and R. Roy, Special functions, Encyclopedia of Mathematics and 
its Applications, 71, Cambridge Univ. Press, Cambridge, 1999. 

[2] T. M. Atanackovic and B. Stankovic, On a numerical scheme for solving differential equa- 
tions of fractional order, Mech. Res. Comm. 35 (2008), no. 7, 429-438. 

[3] P. L. Butzer, A. A. Kilbas and J. J. Trujillo, Mellin transform analysis and integration by 
parts for Hadamard-type fractional integrals, J. Math. Anal. Appl. 270 (2002), no. 1, 1-15. 

[4] P. L. Butzer, A. A. Kilbas and J. J. Trujillo, Stirling functions of the second kind in the 
setting of difference and fractional calculus, Numer. Funct. Anal. Optim. 24 (2003), no. 7-8, 
673-711. 

[5] J. Hadamard, Essai sur l'etude des fonctions donnees par leur developpement de Taylor, 
Journ. de Math. 4 (1892), no. 8 , 101-186. 

[6] U. N. Katugampola, New approach to a generalized fractional integral, Appl. Math. Corn- 
put. 218 (2011), no. 3, 860-865. 

[7] A. A. Kilbas, Hadamard-type fractional calculus, J. Korean Math. Soc. 38 (2001), no. 6, 
1191-1204. 

[8] A. A. Kilbas, H. M. Srivastava and J. J. Trujillo, Theory and applications of fractional 
differential equations, North-Holland Mathematics Studies, 204, Elsevier, Amsterdam, 2006. 



15 



[9] A. A. Kilbas and A. A. Titioura, Nonlinear differential equations with Marchaud- 
Hadamard-type fractional derivative in the weighted space of summable functions, Math. 
Model. Anal. 12 (2007), no. 3, 343-356. 

[10] K. S. Miller and B. Ross, An introduction to the fractional calculus and fractional differential 
equations, A Wiley-Interscience Publication, Wiley, New York, 1993. 

[11] I. Podlubny, Fractional differential equations, Mathematics in Science and Engineering, 198, 
Academic Press, San Diego, CA, 1999. 

[12] D. Qian, Z. Gong and C. Li, A generalized Gronwall inequality and its application to 
fractional differential equations with Hadamard derivatives, 3rd Conference on Nonlinear 
Science and Complexity (NSC10), Cankaya University, 1-4, Ankara, Turkey, 28-31 July, 
2010. 

[13] S. G. Samko, A. A. Kilbas and O. I. Marichev, Fractional integrals and derivatives, trans- 
lated from the 1987 Russian original, Gordon and Breach, Yverdon, 1993. 



16 



